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Abstract. We calculate diffractive dijet production in deep-inelastic scattering at 
next-to-leading order of perturbative QCD, including contributions from direct and 
resolved photons, and compare our predictions to preliminary data from the HI 
collaboration at HERA. We study how the cross section depends on the factorization 
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breaking. The strong M 7 -dependence, which is present when only the resolved cross 
section is suppressed, is tamed by introducing the suppression also into the initial-state 
NLO correction of the direct part. 
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1. Introduction 

It is well known that in high-energy deep-inelastic ep-collisions a large fraction of the 
observed events are diffractive. These events are defined experimentally by the pres- 
ence of a forward-going system Y with four-momentum py, low mass My (in most 
cases a single proton or a proton plus low-lying nucleon resonances), small momentum 
transfer squared t = (p — py) 2 , and small longitudinal momentum transfer fraction 
X P = q{p ~ Py)/<1P from the incoming proton with four-momentum p to the system 
X (see Fig. |TJ). The presence of a hard scale, as for example the photon virtuality 
Q 2 = —q 2 in deep-inelastic scattering (DIS) or the large transverse jet momentum p^ in 
the photon-proton center-of-momentum frame, should then allow for calculations of the 
production cross section for the central system X with the known methods of perturba- 
tive QCD. Under this assumption, the cross section for the inclusive production of two 
jets, e + p — > e + 2jets + X' + Y, can be calculated from the well-known formulae for 
jet production in non- diffractive ep collisions, where in the convolution of the partonic 
cross section with the parton distribution functions (PDFs) of the proton the latter ones 
are replaced by the diffractive PDFs. In the simplest approximation, they are described 
by the exchange of a single, factorizable pomeron/Regge-pole. 

The diffractive PDFs have been determined by the HI collaboration at HERA from 
high-precision inclusive measurements of the DIS process ep — *> eXY using the usual 
DGLAP evolution equations in leading order (LO) and next-to-leading order (NLO) 
and the well-known formula for the inclusive cross section as a convolution of the in- 
clusive parton- level cross section with the diffractive PDFs pQ. For a similar analysis 
of the inclusive measurements of the ZEUS collaboration see |2j. For inclusive diffrac- 
tive DIS it has been proven by Collins that the formula referred to above is applicable 
without additional corrections and that the inclusive jet production cross section for 
large Q 2 can be calculated in terms of the same diffractive PDFs |3j. The proof of this 
factorization formula, usually referred to as the validity of QCD factorization in hard 
diffraction, also appears to be valid for the direct part of photoproduction (Q 2 ~ 0) 
or low-Q 2 electroproduction of jets jSj- However, factorization does not hold for hard 
processes in diffractive hadron-hadron scattering. The problem is that soft interactions 
between the ingoing two hadrons and their remnants occur in both the initial and final 
state. This agrees with experimental measurements at the Tevatron Predictions 
of diffractive dijet cross sections for pp collisions as measured by CDF using the same 
PDFs as determined by HI pQ overestimate the measured cross section by up to an 
order of magnitude 4j. This suppression of the CDF cross section can be explained 
by considering the rescattering of the two incoming hadron beams which, by creating 
additional hadrons, destroy the rapidity gap 0. 

Processes with real photons (Q 2 ~ 0) or virtual photons with fixed, but low Q 2 
involve direct interactions of the photon with quarks from the proton as well as re- 
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Figure 1. Diffractive scattering process ep — » eXY, where the hadronic systems X 
and Y are separated by the largest rapidity gap in the final state. 

solved photon contributions, leading to parton-parton interactions and an additional 
remnant jet coming from the photon (for a review see [Hj). As already said, factoriza- 
tion should be valid for direct interactions as in the case of DIS, whereas it is expected 
to fail for the resolved process similar as in the hadron-hadron scattering process. In 
a two-channel eikonal model similar to the one used to calculate the suppression factor 
in hadron-hadron processes [Sj, introducing vector-meson dominated photon fluctua- 
tions, a suppression by about a factor of three for resolved photoproduction at HERA 
is predicted jjj. Such a suppression factor has recently been applied to diffractive di- 
jet photoproduction [HI El and compared to preliminary data from HI [TO] and ZEUS 
[TT] . While at LO no suppression of the resolved contribution seemed to be necessary, 
the NLO corrections increase the cross section significantly, showing that factorization 
breaking occurs at this order for resolved photoproduction and that a suppression factor 
R = 0.34, in agreement with the prediction of [7j, gives a reasonable description of the 
experimental data of ^0] and [TT] . 

As already mentioned in our earlier work [HI E] , describing the factorization break- 
ing in hard photoproduction as well as in electroproduction at very low Q 2 ^2] by 
suppressing the resolved contribution only may be problematic. An indication for this 
is the fact that the separation between the direct and the resolved process is uniquely 
defined only in LO. In NLO these two processes are related. The separation depends on 
the factorization scheme and the factorization scale My. The sum of both cross sections 
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is the only physically relevant cross section, which is approximately independent of the 
factorization scheme and scale [T^]. We demonstrated already in (HJ |H| that by multi- 
plying the resolved cross section with the suppression factor R = 0.34, the correlation 
of the M 7 -dependence between the direct and resolved part is destroyed and the sum 
of both parts has a stronger M 7 -dependence than for the unsuppressed case (R — 1), 
where the M 7 -dependence of the NLO direct cross section is compensated to a high 
degree against the M 7 -dependence of the LO resolved part. The introduction of the re- 
solved cross section is dictated by perturbation theory. At NLO, collinear singularities 
arise from the photon initial state, which are absorbed at the factorization scale into 
the photon PDFs. This way the photon PDFs become M 7 -dependent. The equivalent 
M 7 -dependence, just with the opposite sign, is left in the NLO corrections to the direct 
contribution. With this knowledge, it is obvious that we can obtain a physical cross 
section at NLO, i.e. the superposition of the NLO direct and LO resolved cross section, 
with a suppression factor R < 1 and no M 7 -dependence left, if we also multiply the 
In M 7 -dependent term of the NLO correction to the direct contribution with the same 
suppression factor as the resolved cross section. 

It is the purpose of this paper to present the special form of the In M 7 -term in the 
NLO direct contribution and then, after multiplying it with the suppression factor R, 
to demonstrate that the M 7 -dependence of the physical cross section cancels to a large 
extent in the same way as in the unsuppressed case (R— 1), if we modify the suppres- 
sion mechanism in the way as explained above. Of course, some small M 7 -dependence 
will remain, also in the case with the suppression of the In M 7 -part in the direct contri- 
bution, which is due to the evolution of the photon PDFs and/or the NLO correction 
to the resolved part, which is left uncompensated as long as the NNLO contribution to 
the direct part is not included. 

In addition to checking the M 7 scale dependence we shall study how the NLO cross 
section depends on the factorization scheme for the photon PDFs. These studies can 
be done for photoproduction (Q 2 ~ 0) as well as for electroproduction with fixed, small 
Q 2 . Since in electroproduction the initial-state singularity in the limit Q 2 — > is more 
directly apparent than for the photoproduction case, where this singularity appears as a 
pole in d — 4 using dimensional regularization, we shall consider in this paper the low-Q 2 
electroproduction case just for demonstration. This diffractive dijet cross section has 
been calculated recently [T2j. In this work da/dQ 2 as a function of Q 2 and da/dzp as 
a function of zjp for various Q 2 + p^ bins was calculated for HI kinematical conditions 
and compared to the preliminary experimental data from HI [Tl] . Here, zp is the par- 
ton momentum fraction in the pomeron. In the following we shall consider these cross 
sections with the same kinematical constraints as in [T2] . 

A consistent factorization scheme for low-Q 2 virtual photoproduction has been de- 
fined and the full (direct and resolved) NLO corrections for inclusive dijet production 
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have been calculated in jTHj. In this work we adapt this inclusive NLO calculational 
framework to diffractive dijet production at low-Q 2 in the same way as in [12] . except 
that we multiply the In M 7 dependent terms as well as the resolved contributions with 
the same suppression factor R = 0.34 as in our earlier work [SI El El- The exact value of 
this suppression factor may change in the future, when better data for photoproduction 
and \ow-Q 2 electroproduction have been analyzed. 

In the next Section, we shall define the HI kinematical constraints and present the 
In M 7 -dependent terms in the NLO corrections to the direct cross section, as they are 
reported in ^H]- We shall also define which specific contribution of this initial-state NLO 
correction will be suppressed with the same factor R as the resolved part. Then we give 
our results concerning the scheme dependence of the photon PDFs by looking at da/dQ 2 
and da/dzjp using two different schemes. After this we present the In M 7 - dependence of 
the partly suppressed NLO direct and the fully suppressed NLO resolved cross section 
da/dQ 2 and their sum for the lowest Q 2 bin. After this we study the M 7 -dependence 
of da/dQ 2 for the other Q 2 -bins and of da/dzjp and compare always with the HI data. 
As the last point we study the contribution of the hadron-like or vector-dominance part 
of the resolved cross section in order to see whether the suppression of this part alone 
could be an alternative to solve the M 7 scale problem as suggested in [SJ E] • 

2. Kinematical constraints and results 

The factorization scheme for virtual photoproduction has been defined and the full NLO 
corrections for inclusive dijet production have been calculated in They have been 
implemented in the NLO Monte Carlo program JETVIP ^E]- We have adapted this 
NLO framework to diffractive dijet production. According to JTHj, the subtraction term, 
which is absorbed into the PDFs of the virtual photon f a / 7 (x, M 7 ), is of the form 



Qi is the fractional charge of the quark qi, p\ and P2 are the momenta of the two 
outgoing jets, and po and q are the momenta of the ingoing parton and virtual photon, 
respectively. Since Q 2 <C M 2 , the subtraction term in Eq. (0) is large and is therefore 
resummed by the DGLAP evolution equations for the virtual photon PDFs. After the 
subtraction of Eq. (fl|). the finite term, which remains in the matrix element for the NLO 
correction to the direct process, is [THj 



This term is defined by the requirement that it is equal to for real photons, if 

Q 2 = 0. For our study of the M 7 -dependence of the physical cross section it is essential to 




(1) 



(2) 




(3) 
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recognize that the M 7 -dependence in Eq. (J2J) is the same as in the subtraction term Eq. 
(JIJ, i.e. lnM 7 is multiplied with the same factor (except that the subtraction term S has 
to be multiplied by — l/(2iV c )). As already mentioned, this yields the M 7 -dependence 
before the evolution is turned on. In the usual non-diffractive dijet photoproduction 
these two M 7 -dependences, i.e. in Eqs. ((TJ) and (jHJ), cancel, when the NLO correction to 
the direct part is added to the LO resolved cross section [T3|. Then it is obvious that the 
approximate M 7 -independence is destroyed, if the resolved cross section is multiplied by 
the suppression factor R to account for the factorization breaking in the experimental 
data. To remedy this deficiency, we propose to multiply the In M 7 -dependent term in 
M(<5 2 )ms w ith the same suppression factor as the resolved cross section. This is done 
in the following way: We split M(Q 2 )j^ into two terms using the scale p^ in such 
a way that the term containing the slicing parameter y s , which was used to separate 
the initial-state singular contribution, remains unsuppressed. In particular, we replace 



where R is the suppression factor. This expression coincides with M{Q 2 )y^ in Eq. (J3j) 
for R = 1, as it should, and leaves the second term in Eq. (J1J) unsuppressed. In Eq. (J2J 
we have suppressed in addition to ln(M 7 /p^ 2 ) also the z-dependent term hi(z/ (1 — z)), 
which is specific to the MS subtraction scheme as defined in ^H]- To keep this term 
unsuppressed, i.e. to move it to the second term in Eq. (jlj), would again produce an in- 
consistency with the suppressed resolved contribution. The second term in Eq. (j3J) must 
be left in its original form, i. e. being unsuppressed, in order to achieve the cancellation 
of the slicing parameter dependence of the complete NLO correction in the limit of very 
small Q 2 or equivalently very large s. With these arguments the splitting of the terms 
in Eq. ©, which can be suppressed, and those terms, which remain unsuppressed, is 
almost unique. It is clear that the suppression of this part of the NLO correction to 
the direct cross section will change the full cross section only very little as long as we 
choose M 7 ~ p* T . The first term in Eq. (0J), which has the suppression factor R, will be 
denoted by DIRis in the following. 

To study the left-over M 7 -dependence of the physical cross section, we have 
calculated the diffractive dijet cross section with the following kinematical constraints as 
in the HI experiment [T3|: The electron and proton beam energies are 27.5 and 820 GeV, 
respectively, 4 < Q 2 < 80 GeV 2 is the range of the squared electron momentum transfer, 
0.1 < y < 0.7 is the interval of virtual photon energy fraction, xp < 0.05, \t\ < 1 GeV 2 , 
My < 1.6 GeV, PT.jeti,2 > 5(4) GeV, and jet rapidities —3 < f]jeti,2 < 0- J e ^ s are defined 
by the CDF cone algorithm with jet radius equal to one. The asymmetric cuts for the 
transverse momenta of the two jets are required for infrared stable comparisons with the 
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NLO calculations ^7]. The original HI analysis actually used a symmetric cut of 4 GeV 
on the transverse momenta of both jets ^H]- The data have, however, been reanalyzed 
for asymmetric cuts [H|. In the extraction of the diffractive parton densities the usual 
Regge factorization ansatz [TH] 



has been employed, which we shall use also in our computations. This ansatz neglects 
a possible scale dependence of fjp/ p , which could be responsible for a change of the xjp 
dependence between the soft and the hard pomeron exchange. In the analysis of HI 
the intercept ctjp(O) is extracted from the diffractive DIS data. For an alternative 
approach, where the factorization assumption in Eq. (JHJ) is modified see [20J. 

For the NLO resolved virtual photon predictions, we have used the PDFs SaSID 
|2~T] and transformed them from the DIS 7 to the MS scheme as in [TS]. If not stated 
otherwise, the renormalization and factorization scales at the pomeron and the photon 
vertex are equal and fixed to p* T = Pr,jeti- We include four flavours, i.e. nf = 4 in the 
formula for a s and in the PDFs of the pomeron and the photon. The SaS-type PDFs 
for the virtual photon are chosen, since in their parametrization the scale dependence 
following from the evolution is given separately for the hadronic part and the point-like 
part, which is called the anomalous part in [2~T] . 

With these assumptions we have calculated the same cross section as in our pre- 
vious work [12]. First we investigated how the cross section da/dQ 2 depends on the 
factorization scheme of the PDFs for the virtual photon, i.e. da/dQ 2 is calculated for 
the choice SaSID and SaSlM. Here da/dQ 2 is the full cross section (sum of direct and 
resolved) integrated over the momentum and rapidity ranges referred to above. The re- 
sults are shown in Fig. where da/dQ 2 with SaSID (full line) is compared to da/dQ 2 
with SaSlM (dashed line). The theoretical cross sections are integrated over the same 
Q 2 bins as the experimental ones. In this figure and in all the following ones, except 
Fig. HI where we study the explicit M 7 -dependence, both predictions are obtained with 
suppressed resolved cross section and the additional suppression of DIRjs as described 
above with suppression factor R = 0.34. As we can see in Fig. the choice of the 
factorization scheme of the virtual photon PDFs has negligible influence on da/dQ 2 for 
all considered Q 2 . The predictions agree reasonably well with the preliminary HI data 
|14j . Comparing with our previous results ^2] we notice that the cross sections in Fig. 
12] are somewhat smaller than the corresponding ones in Fig. 2 of jT2|. This is, however, 
not the effect of the additional suppression of a NLO term in the direct cross section. 
The difference comes instead from a different choice of the overall scales in J2] which 
was fixed to the average (p^) = v^40 GeV. Another check on the effect of a different 
choice of the factorization scheme is shown in Fig. El where the dependence of the cross 
section da/dzjp as a function of zp is plotted for various ranges of Q 2 +p^. From this 
plot we see that the scheme dependence is again very small for all Q 2 + p*^ ranges and 




(5) 
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Figure 2. Dependence of the diffractive dijet cross section at HERA on the 
squared photon virtuality Q 2 . Preliminary HI data |14) are compared with theoretical 
NLO predictions including resolved virtual photon and direct contributions with a 
suppression factor of R = 0.34 for the resolved and for the DIRis term using virtual 
photon PDFs SaSID (full) and SaSIM (dashed) from [23] . 

for all zp G [0; 1]. Except for the lowest zjp bin and for the two highest Zjp bins in 
the two upper Q 2 + ranges, we observe reasonable agreement with the data. The 
agreement is best for the smallest Q 2 + G [20; 35] GeV 2 region and Zjp > 0.2. 

We now turn to the M 7 -dependence of the cross section with a suppression factor 
for DIRis, which is the main part of this paper. To show this dependence for the two 
suppression mechanisms, (i) suppression of the resolved cross section only and (ii) ad- 
ditional suppression of the DIRjg term as defined in Eq. (J3J) in the NLO correction of 
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Figure 3. Dependence of the diffractive dijet cross section at HERA on the parton 
momentum fraction in the pomeron zp for different ranges of Q 2 + p^ 2 . Preliminary 
HI data ^3] are compared with two theoretical predictions (as in Fig.EJ) using SaSID 
(full) and SaSIM (dashed) virtual photon PDFs [23]. 



the direct cross section, we consider da/dQ 2 for the lowest Q 2 -bin, Q 2 6 [4,6] GeV 2 . In 
Fig. HI this cross section is plotted as a function of £ = M 7 /p^ in the range £ G [0.25,4] 
for the cases (i) (light full curve) and (ii) (full curve). We see that the cross section 
for case (i) has an appreciable ^-dependence in the considered £ range of the order of 
40%, which is caused by the suppression of the resolved contribution only. With the 
additional suppression of the DIRis term in the direct NLO correction, the ^-dependence 
of da/dQ 2 is reduced to approximately less than 20%, if we compare the maximal and 
the minimal value of da/dQ 2 in the considered £ range. The remaining ^-dependence is 
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Figure 4. Photon factorization scale dependence of resolved and direct contributions 
to da/dQ 2 together with their weighted sums for (i) suppression of the resolved cross 
section and for (ii) additional suppression of DIRis, using SaSID virtual photon PDFs 

EH- 



caused by the NLO corrections to the suppressed resolved cross section and the evolution 
of the virtual photon PDFs. How the compensation of the M 7 -dependence between the 
suppressed resolved contribution and the suppressed direct NLO term works in detail 
is exhibited by the dotted and dashed-dotted curves in Fig. |U The suppressed resolved 
term increases and the suppressed direct NLO term decreases by approximately the 
same amount with increasing £. In addition we show also da/dQ 2 in the DIS theory, 
i.e. without subtraction of any InQ 2 terms (dashed line). Of course, this cross section 
must be independent of £. This prediction agrees very well with the experimental point, 
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Figure 5. Same as in Fig. for M 7 = (full)and M 1 = kp* T (dashed) and 

comparison with preliminary HI data using SaSID virtual photon PDFs |21) . 



whereas the result for the subtracted and suppressed theory (full curve) lies slightly 
below. We notice, that for M 7 = p* T the additional suppression of DIRis has only a 
small effect. It increases da/dQ 2 by 5% only. 

In order to get an idea about the M 7 scale dependence of da/dQ 2 for the other Q 2 
bins and for dcr /dzjp in the Q 2 +p^ ranges as in Fig. El we have computed these two cross 
sections for two choices of M 7 , namely M 7 = p^/4 and M 7 = 4p^ corresponding to the 
lowest and highest £ in Fig.HJ The result for the da/dQ 2 is shown in Fig.El We see that 
the M 7 - dependence in the considered range decreases with increasing Q 2 . This is to be 
expected since the resolved contribution diminishes with increasing Q 2 , so that the NLO 
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Figure 6. Same as in Fig. 03 for M 7 — (full) and M 7 = \p* T (dashed) and 

comparison with preliminary HI data using SaSID virtual photon PDFs |21) . 



corrections to the resolved cross section and the effect of the evolution of the photon 
PDF diminish as well. Similar conclusions can be drawn from Fig. El where we plotted 
da/dzjp in the four Q 2 + py 2 -ranges. For all four ranges the dependence on M 7 is small 
if zjp > 0.4. Only for the two lowest Q 2 + p5?-ranges this cross section depends on M 7 
for the two lowest zp-bins and is strongest for zjp G [0; 0.2] and Q 2 +p^ G [20; 35] GeV 2 . 

As is well known the photon PDFs consist of a point-like or anomalous part, the 
hadron-like part and the gluon part. The compensation of the M 7 -dependence between 
the NLO direct and the LO resolved cross section occurs via the anomalous part. This 
means that this part of the PDFs is closely related to the direct cross section. If one 
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assumes that the direct part obeys factorization , i.e. has no suppression factor, it 
was suggested in 0, that one possibility to avoid the non-compensation of the ./In- 
dependence between the suppressed resolved cross section and the unsuppressed direct 
cross section would be to suppress only the hadron-like and the gluon part of the PDFs. 
Fortunately, in the SaS-type photon PDFs the hadron-like and the anomalous part and 
the effect of their evolution are presented separately, so that the effect of the suppression 
of the hadron-like part can be investigated for all scales. It turns out, however, that the 
hadron-like part contributes only a very small fraction to the total resolved cross section. 
For Q 2 = 5 GeV 2 , this fraction is 3 • 10~ 3 , and it decreases strongly with increasing Q 2 . 
This means that the anomalous component dominates the total resolved cross section 
and a suppression of the hadron-like part alone would not be sufficient to account for 
the experimental data at low Q 2 . Although in the case of photoproduction (Q 2 ~ 0) 
the hadron-like part is supposed to be larger, it would be quite artificial to have two 
different mechanisms for the M 7 -scale compensation depending on whether Q 2 > or 
Q 2 ~ 0. For this reason we prefer the mechanism with the suppression in the DIRis 
contribution to the NLO corrections of the direct cross section, which works for all Q 2 . 

3. Conclusions 

In summary, we propose in this paper a new factorization scheme for diffractive pro- 
duction of jets in low-Q 2 deep inelastic scattering. By suppressing not only the resolved 
photon contribution, but also the unresummed logarithm as well as scheme- dependent 
finite terms in the NLO direct initial state correction, factorization scheme and scale 
invariance is restored up to higher order effects, while at the same time the cut-off in- 
variance required in phase space slicing methods is preserved. 

For pedagogical reasons, we have chosen in this paper the kinematical region of 
finite, but low photon virtuality Q 2 , which exposes and regularizes a logarithmic vir- 
tual photon initial state singularity. We do, however, not rely on the finiteness of 
Q 2 , but rather separate suppressed and unsuppressed terms using the hard transverse 
momentum scale pj,, so that our scheme is equally valid for real photoproduction. Fur- 
thermore, due to the universality of the initial state singularity, our new factorization 
scheme should apply to any other diffractive photon-induced process with a second hard 
scale, such as inclusive hadron-production at large ptp or heavy quark and quarkonium 
production. 

The scheme- and scale invariance has been demonstrated numerically using the 
kinematics of a recent HI analysis, differential in Q 2 and parton momentum fraction in 
the pomeron zjp. Very good stability with respect to scheme- and scale variations and 
good agreement with the experimental data has been found. 

Individual studies of the anomalous and hadronic component of the resolved 



Factorization scheme and scale dependence in diffractive dijet production 



14 



virtual photon show that the parton densities in the virtual photon are dominated 
by the resummed pointlike component, as expected from non- diffractive deep inelastic 
scattering. While suppressing only the hadronic component in the photon PDFs 
would also be scheme- and scale- invariant and therefore be theoretically consistent, 
our numerical study shows that this alternative is thus phenomenologically not viable. 
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